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Abstract. A positive unit point charge approaching from infinity a perfectly spherical isolated 
conductor carrying a total charge of +1 will eventually cause a negatively charged spherical 
cap to appear. The determination of the smallest distance p(d) (d is the dimension of the unit 
sphere) from the point charge to the sphere where still all of the sphere is positively charged is 
known as Gonchar's problem. Using classical potential theory for the harmonic case, we show 
that 1 + p(d) is equal to the largest positive zero of a certain sequence of monic polynomials of 
degree 2d — 1 with integer coefficients which we call Gonchar polynomials. Rather surprisingly, 
p(2) is the Golden ratio and p(4) the lesser known Plastic number. But Gonchar polynomials 
have other interesting properties. We discuss their factorizations, investigate their zeros and 
present some challenging conjectures. 



1. Introduction 

Let S d denote the unit sphere in the Euclidean space Suppose that it is insulated and 

has a total positive charge of +1. In the absence of an external field the charge will distribute 
uniformly with respect to the normalized surface area measure (unit Lebesgue measure) ad- Now 
we introduce a positive unit point charge exterior to the sphere that repels the charge on the 
sphere in accordance with the Newton potential l/r d_1 , where r represents the distance between 
point charges. If this point charge is very close to the sphere, then one would expect it to cause 
a negatively charged spherical cap to appear, while if the point charge is very far from the sphere 
its influence is negligible and the charge on the sphere will be everywhere positive and nearly 
uniformly distributed over the entire sphere. We consider the following question: what is the 
smallest distance from the unit point charge to S d such that the distribution of the positive charge 
on the sphere covers all of the sphere? We will denote this critical distance by p(d). As we shall 
show, l+p(d) equals the largest positive root of the following polynomial equation of degree 2d — 1: 

(1) G(d;z):=\(z-l) d - z-llz*- 1 + (z-lf = 0, d= 1,2,3,.... 

As the question above was communicated to the authors by A. A. Gonchar, we shall refer to 
G(d; z) as Gonchar polynomials. 

Rather surprisingly, p{2) turns out to be the Golden ratio, which is the limit of the ratio of 
successive terms in the Fibonacci sequence F n — F n -\ + F n -2, Fo = Fi = 1 & n d p(4) is the 
so-called Plastic number [10], which is the limit of ratios of successive terms for the less known 
Padovan sequence P n = P n -2 + P n -3, P b = Pi = Pi = 1 (sequence A000931 in Sloane's OEIS [9]). 
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In addition to this curious coincidence, the polynomials z) exhibit rather fascinating prop- 
erties with regard to their irreducibility over the ring of polynomials with integer coefficients, as 
well as the asymptotic behavior of their zeros. Our goal in the next section is to show how the 
Gonchar polynomials are derived and then, in Section 3, to explore some of their properties and 
draw the reader's attention to some related conjectures. 

2. Signed Equilibrium 

A general charge distribution on §> d = {x £ R d+1 : |x| = 1} will be modeled by a signed 
measure n supported on S d with r/(S d ) = 1. The corresponding Newtonian potential and energy 
are given by 

V"(x) = J |x _ 1 y|d _ 1 ^(y), Z(„) = J t™(x)d„(x). 

For example, V(x) is the familiar Coulomb potential when d — 2. For n = aa we expect from the 
rotational invariance that the potential V' Td is constant on S d and on concentric spheres. This is 
a well-known fact in potential theory, but for completeness we provide a proof below. 

Lemma 1. The potential V ad (x) satisfies the following relations 

'\ for |x|<l, 

l/lx^- 1 for |x| > 1. 



(2) V ad (x) 



Proof. Indeed, for any fixed x the Newtonian kernel /(y) = |x — yl 1 d is a harmonic function 
in R d+1 \ {x}, that is A/(y) = J2i=i fym (y) = 0- Therefore, by the mean value property for 
harmonic functions we have 

(3) V™(x) = / /(y) ^(y) = /(0) = r 4rr for |x| > 1. 



Considering the sequence of functions f n {y) = + l)x/n — y| 1_d and applying the monotone 
convergence theorem allows us to extend this to the unit sphere giving V r<Td (x) = 1 for |x| = 1. 
By the maximum principle for harmonic functions, the potential V ad is constant everywhere in 
the closed unit ball ("Faraday cage effect"), since it assumes its extreme values on the boundary 
§ d , where it is constant. (One can also verify this directly by applying the identity |x — y| 2 = 
|x| 2 |^-y| 2 and(3).) □ 

In the classical Coulomb case (d = 2), a standard electrostatics problem (see [3, Ch. 2]) is 
to find the charge distribution on a charged, insulated, conducting sphere in the presence of an 
external field (such as generated by a positive point charge q off the sphere). This motivates the 
following definition. 

Definition 2. Given a continuous function Q(x) on § d , we call a signed measure ijq supported on 
S d and of total charge VQ(§> d ) = 1 a signed equilibrium associated with Q if its weighted Newtonian 
potential is constant on § d ; that is, for some constant Fq, 

V m (x) + Q(x) = F Q everywhere on S d . 

In the above definition, the function Q is referred to as an external field. Thus Lemma 1 
establishes that o~d is the (signed) equilibrium measure on S d in the absence of an external field. The 
uniqueness of the signed equilibrium is settled by the following known proposition [1, Lemma 23], 
whose proof is presented in the Appendix. 

Proposition 3. If a signed equilibrium t]q associated with an external field Q on E> d exists, then 
it is unique. 
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Here, we are concerned with the external field generated by a positive point charge of amount 
q located at a — Rp, R > 1, where p is the North Pole of S d ; that is, p = (0, 0, . . . , 0, 1). Such a 
field is given by 

(4) q„ (x):= __|_ xer m 

and the associated signed equilibrium is described in the next result. (For the general result, see 
WO 

Lemma 4. The signed equilibrium t)q associated with the external field Q = Qa.,q of (4) is abso- 
lutely continuous with respect to the (normalized) surface area measure on E d ; that is, dr]Q(x) = 
r]' R q (x)dad(x), and its density is given by 

For the classical Coulomb case this relation is well-known from elementary physics (cf. [3, 
p. 61]). 

Proof. We note first that if |x| > 1, then the function u(z):=|z — x| 1_d is harmonic in |z| < 1, and 
since the Poisson integral formula [2] preserves harmonic functions, we have 



|z|<l. 



1 f 1 1 - Izl 2 

( 6 ) i id^T = / i id=Ti id+T dcTd (y)> 

|z-x| d 1 7g d |y-x| d 1 |z — y| d+1 

Using a monotone convergence theorem argument as in Lemma 1 we extend (6) to |x| = 1. We 
shall make use of the identity \Ry — x| = |y — Rx\ for any |x| = |y| = 1. Then from (6) we obtain 
with z = p/R and |x| = 1, 
(7) 

n( x = g = q = g = f q i - |z| 2 . . 

VW l^p-x^- 1 \p-Rx\d- 1 ^-^z-xl^- 1 J sd R d - 1 \y-x\ d - 1 \z-y\ d + 1 d[y > 



|y -x\^\p-Ry\™ dadiy) = I ly-x^ |a-y|^ ^ (y) ' 
The Poisson integral formula applied to the constant function u(z) = qjR! 1-1 and z = p/R yields 

q f q(R 2 - 1) 

This together with (7) implies that for dv:=r\ R q dod we have 



/ 



du = l and ^(x) + Q(x) = l + -| ZT , x e § d , 



which proves the lemma. □ 

Remark. A surprising aspect of the equilibrium support, which varies with R, is illustrated in 
Figure 1 for § 2 and q = 1. Writing r]Q a l = VQ = Vq ^ Vq i n terms of its Jordan decomposition, it 
is clear on geometrical and physical grounds that the support of rfq is a spherical cap centered at 
the South Pole, while its complement (a spherical cap centered at the North Pole) is the support 
of t]q, see Figure 1(a). Although somewhat counter-intuitive, we see from Figure 1(b) that, as 
the unit charge approaches the sphere (R — > 1), from a certain distance on, the support of 7]q 
occupies increasingly more of the sphere. This phenomenon can be explained by the geometry of 
the sphere, where an increasingly 'needle-like' negative sink forming at the North Pole is balanced 
by an increasingly uniform positive part so that the total charge is always one. 

From formula (5) we observe that the minimum value of the density rj' R q (x) is attained at the 
North Pole p, and its value there is 

(8) W q q ( R2 ~^ -1 | « *( R2 - 1 ) 



R*- 1 |a-p| d+1 R*- 1 (R-l) d+1 ' 
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This leads to the following theorem. 

Theorem 5. For the external field Q a ,<j( x ) of (4) with a = Rp, the signed equilibrium is a positive 
measure on all ofS d if and only if R> R q , where R q is the unique (real) zero in (1, +oo) of the 
polynomial 



(9) 



G(d,q;z):-- 



l) d /q- 



1 



- 1 + (z-l)° 



In particular, the solution to Gonchar's problem is given by p(d) = R\ — 1. 

Proof. We use the fact that the support of r\Q & q is all of S d if and only if n' R (p) > in (8) or 
equivalently that G(d,q;R) > 0. Hence we seek the number R > 1 such that G(d,q;R) = 0. 
Observing that G(d, q; 1) < and G(d, q;x) > for x > 1 sufficiently large, there exists at least 
one value R q such that G{d, q; R q ) = 0. Moreover, this root in (1, oo) is unique as can be seen by 
applying Descartes' Rule of Signs to 

d-l / > ,\ d— 1 



G(d,q;l 



- w 



q ^ 



d-l 

m 



E 

m=0 



d-l 

m 



w 



□ 



Curiously, for d = 2 and c? = 4 we obtain the Golden ratio and the Plastic number as answers 
to Gonchar's problem as mentioned in the Introduction. Furthermore, for large values of d the 
asymptotic analysis (provided in the Appendix) shows that 



(10) 



Ra 



[log(3 g )]/d + 0(l/d 2 ) 



as d 



oo. 



The appearance of 2 as the limit as d — > oo of the critical distances R q (q > fixed) can be 
explained by studying the asymptotic behavior of G(d,q;R) as d — > oo. Indeed, for R > 2 + e 
(e > fixed) it goes to +oo as d — > oo, while for 1 < i? < 2 — e(e>0 fixed) it goes to — oo, 
leaving the number 2 as the only candidate. 



3. The polynomials G(d;z) 

In the following we investigate the family of Gonchar polynomials G(d, z) = G(d, 1; z) given in 
(1). Aside from the solution to Gonchar's problem, these polynomials are interesting in themselves 
and their distinctive properties merit further study. 
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A polynomial P with real coefficients is called (self-)reciprocal if its reciprocal polynomial 
P* (z):~z dcg p P(l I z) coincides with P{z). In other words, the coefficients of z k and of z dcgP ~ k 
in P(z) are the same. Notice that G(d; z) is self-reciprocal for even d since 



z 2d - 1 G{d;l/z)= (l-z) d -z-l z d - l + {l-z 

Consequently, if C, is a zero of G(d; z), then so is 1/Q for even d. For odd d we infer from 

G(d; z) + G*{d; z) = [l + {-l) d ] (z - i f (z^ 1 + l) - 2 (z d + z^ 1 ) = -2 (z d + z^ 1 ) 

that the coefficients of z k and z 2d ~ 1 ~ k in G(d; z) sum to zero except for the 'innermost' pair. 

3.1. Factorizations and Irreducibility. With the aid of symbolic computation programs one 
can find factorizations and check irreducibility of explicitly given polynomials. For d = 1, 2, . . . , 7, 
we thereby obtain the following factorizations over the integers of G(d; z). 



vd-l i ft „\ d 



G(l;z 
G(2;z 
G(3;z 
G(4;z 
G(5;z 
G(6;z 
G(7;z 



z- 3, 

(z + 1) (z 2 -3z + l) , 

z 5 -3z 4 + 3z 3 -5z 2 + 3z-l, 

{z + 1) (z 3 - 3z 2 +2^-1) (z 3 - 2z 2 + 3z - l) , 

z 9 - 5z 8 + 10z 7 - 10z 6 + 5z 5 - 7z 4 + 10z 3 - 10z 2 + 5z - 1, 

(z + 1) (z 2 - z + l) (z 8 - 6z 7 + 15z 6 - 21z 5 + 21z 4 - 21z 3 + 15z 2 - 6z + l) , 

z 13 - 7z 12 + 21z n - 35z 10 + 35z 9 - 21z 8 + 7z 7 - 9z 6 + 21z 5 - 35z 4 + 35z 3 - 21z 2 + 7z - 1. 



One can easily verify that (z + 1) divides G(d; z) if and only if d is even. Furthermore, the 
factor z 2 — z + 1 arises in the following cases. 

Proposition 6. The cyclotomic polynomial z 2 — z + 1 divides G(d; z) if and only if 6 divides d. 

Proof. Note that if ( is a zero of z 2 — z + 1, then ( 2 = ( — 1 and ( 3 = —1. Using formula (1), it 
readily follows that G(d; () = whenever d = (mod 6) and G(d; () ^ otherwise. □ 

A general irreducibility result has so far eluded the authors. The Eisenstein criterion and, 
in general, reduction to finite fields, seem not to be effective tools for studying the polynomials 
G(d; z). However, using Mathcmatica, we verified the following conjecture for d up to 500. 

Conjecture 1. Set £(d; z) = 1 for d odd, £(d;z):=z + 1 for d even but not divisible by 6, and 
£(d;z):=(z + l)(z 2 — z + 1) if 6 divides d. Then G(d;z)/£(d;z) is irreducible over the rationals 
except for d — A, 8 and 12. 

Regarding the exceptional cases we record that, in addition to G(4; z) as given above, 

G(8; z) = (z + 1) (z 4 - 3z 3 + 3z 2 - 3z + l) 

x (z 10 - 6z 9 + 16z 8 - 24z 7 + 24z 6 - 21z 5 + 24z 4 - 24z 3 + 16z 2 - 6z + l) , 

G(12; z) = (z + 1) (z 2 - z + 1) (z 6 - 4z 5 + 5z 4 - 3z 3 + 5z 2 - Az + l) 

x (z 14 - 8z 13 + 29z 12 - 62Z 11 + 85z 10 - 77z 9 + 48z 8 

- 33z 7 + 48z 6 - 77z 5 + 85z 4 - 62z 3 + 29z 2 - 8z + l) . 

3.2. Zeros of Gonchar polynomials. Figure 2 illustrates some features of G(d; z) on the real 
line and suggests some of its general properties. Note the qualitatively different behavior for even 
and odd d. 

Depending on the parity of d, the polynomial equation G(d; z) = has either one (d odd) or 
three (d even) real simple roots. More precisely, the following holds. 

Proposition 7. If d is odd, then G(d; z) has precisely one real zero, which is simple and lies in 
the interval (2,3]. If d is even, then G(d;z) has exactly three real zeros: one at —1, one in the 
interval (1/3, 1/2) and one in the interval (2,3); all these zeros are simple. 
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Proof. Note that for d > 2 there holds: 

G(d;-l) = (-2) d [l-(-l) d ], G(d;0) = (-l) d , G(d; 1) = -2, G(d; 2) - 1 - 2 d . 

The assertion of the proposition is trivial for G(l;z) = z — 3. So let d > 2. Since G(d; 2) < 
and G(d; 3) > 0, the polynomial G(d; z) has at least one real zero in the interval (2, 3) and, as we 
observed in the proof of Theorem 5, this is its only zero on [1, oo) and must be simple. For odd d, 
each of the terms [(a: — l) d — x — X\x d ~ x and (x — l) d is negative for x < 1 and hence so is G(d; x); 
thus G(d; x) has no zeros outside (2, 3). For even d, the self-reciprocity of G(d; x) implies that to 
each zero £(d) of G(d;x) in [l,oo) there is a zero = in (0, 1]. Consequently, by the 

first part of the proof, G(d; x) has one and only one zero in (0, 1) and this zero is simple and lies 
in the interval (1/3, 1/2). It remains only to consider the interval (— oo, 0]. Clearly G(d; —1) = 0, 
and by computing G'(d; x) and analyzing its sign (in particular, G'(d; —1) > 0), one can show that 
G(d;x) is strictly increasing on (— oo, — 1). Thus G(d;x) < on (— oo,— 1) and, by self-reprocity, 
G{d)x) > on (-1,0). □ 

Regarding the behavior as d increases of the zeros in (2,3), M. Lamprecht [4] proved the 
following. 

Proposition 8. The zeros £(d) of G(d; z) in the interval (2,3] form a strictly monotonically 
decreasing sequence with limit point 2 (compare with (10) j. 

We now turn to the study of the complex zeros of G(d;z). Observe that by Proposition 7, 
G(d; z) has either d — 1 (if c? is odd) or d — 2 (if d is even) pairs of complex conjugated zeros 
(counting multiplicity). In Figure 3 we have plotted the zeros of G(d; z) for d = 9, 10, 11, and 12 
along with the two unit circles Co, C\ centered respectively at and 1. Notice that these circles 
intersect at the points (1 ± iy3)/2 and, by Proposition 6, these points are zeros of G(d; z) if and 
only if d is a multiple of 6. The zeros of G(d; z) seem to occur roughly into three categories: zeros 
close to Cq (indicated by x), zeros close to C\ (indicated by +) and zeros close to the vertical line 
x = 1/2 (indicated by ♦). The numbers N\, N2, N3 of zeros in each of the categories are listed in 
Table 1 for d = 1, . . . , 12, 42, from which it appears that these numbers are nearly the same. We 
will discuss this further in Section 3.3 (see Conjecture 2). 

Figure 3 suggests that for even d the polynomial G(d; z) may have zeros lying precisely on Co- 
This, in fact, is the case as we now prove. 

For this purpose it is convenient to define Sq^.—I if 6 divides d and 8o\d'-=Q otherwise. 

Proposition 9. If d is even, G{d\ z) has exactly 4( \_{d — 1)/6J + ^6|d) + 1 zeros on the unit circle 
Co; all are simple and satisfy Re 2: < 1/2. Their positions are determined by the solutions of the 
equation g(9) = f(9) given in (11) below. 

If d is odd, there are no zeros ofG(d;z) on the unit circle Co- 
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Figure 3 and Table 1. Counting the zeros of G(d;z), n = 2d — 1. Two 
additional zeros appear when 6 divides d, see shaded entries in table. 



Proof. Since z = 1 is not a zero of G(d; z), the equation G(d; z) — is equivalent to 

z*- 1 + 1 = 



z-l(z- l) 11 - 1 ' 

Substituting z = e 1 ® and changing to trigonometric functions we arrive at 

(ii) 9 (8y.=(-i) d/2 ^ ( d —±e) = (-i) d/ \" os{0 ! 2 L =-.f(e) 



2 



[2i sin(0/2) 



Suppose d is even. Then f(9) = cos(#/2)/[2sin(#/2)] d is real-valued. Since the complex zeros 
of G{d; z) occur in complex conjugate pairs, we may assume < 8 < ir. Using elementary 
calculus one shows that the function f(9) is monotone decreasing and convex on (0,7r). Hence, 
f(0) > /(7r/3)+/'(7r/3)(6»-7r/3) on(0,7r/3). This implies that f{9) > 1 for 9 e (0, w/3-a), where 
a = (4 — 2v3)/(3d + 1). On the other hand, on [tt/3 — a, 7r/3) one can show that \g(0)\ < f(0). 
Hence, in (11) we have only to consider the range [n/3, n] whereupon f(9) < f{ir/2>) = V3/2 < 1. 

The function g is a cosine function with period Aw/(d — 1), where the sign factor ensures that 
g has a positive derivative at 9 — 7r (which is also an intersection point of / and g). The other 
intersection points of / and g occur in the half-periods where g > 0. On such a half-period / the 
convex function / and the concave function g (when restricted to /) can intersect in at most two 
points counting multiplicity as can be seen by applying Rolle's theorem to the strictly convex 
function f — g. That there are at least two such points on / can be seen from the intermediate 
value theorem applied to the same function (/ — g > at the endpoints of I and / — g < at the 
midpoint of J) as illustrated in Figure 4 for the three canonical cases. 

Since there are |_(cf — 1) /6J full periods oft; in [7r/3,7r] plus one more period (partially contained 
in [7r/3,7r]) whenever 6 divides d, we have 2(\_(d — 1)/6J + <5g|d) zeros in the upper half-plane and 
that many conjugate zeros of G(d; z) in the lower half-plane. 

For odd d, equation (11) has no real solution in (0, 7r]. Since 1 is not a zero of G(d;z), there 
are no zeros of G(d; z) on Co- □ 

Theorem 10. All zeros of G{d; z) are simple for each d > 1. 
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Figure 4. Graphs of f(8) and g(0); typical cases when solving equation (11). 

Proof. The cases d = 1 and d = 2 are obvious, so assume d > 3. Let £ be a zero of G(d; z). Then 
£ is simple if G'(d; £) ^ 0. By means of some helpful substitutions for the expressions in braces*, 
we find that 

C (C - i) G\d- = d( {(c - i) d (C"- 1 + 1) } - C l (C - i) + (d i) { [(C - i) d - c - 1] c"- 1 } (C - i) 
= (d - 1) [c d+1 - (C - l) d+1 ] + (d + 1) c d 
= (d-i)c{c d ~(C-i) d } + (d-i)(C-i) d + (d + i)C d 

= (d - 1) (C - l) d (C d + 1) + 2C d =:P(d; 0- 

Suppose to the contrary that P(d; Q = 0. Then on replacing (C - l) d by -2( d /{(d - l)(( d + 1)] 
in the formula (1) for G(d; £) we get 



(12) 



(C d + 1) G(d; C) = C^ 1 (c d+1 + ^\( d + + l) 



The polynomial Q(d; z) of degree d + 1 obtained by replacing £ by z in the second parenthetical 
expression in (12) has three real zeros (at —1) and d— 2 complex zeros for even d and two negative 
zeros ('near' —1) and d — 1 complex zeros for odd dJ This can be seen from the facts that 



Q(d;-l) = - 2[ \ { y d] , Q'(d;-l)=[l-(-r/] d " ' 



1 



1' 



Q"(d;-1) = 



and Descartes' Rule of Signs. Substituting w = e 1 ^ and using trigonometric functions, we arrive 
at 

Q(d; e 1 ^) = ——. — : cos ^ cos ^ ( d + tan ^ tan ^ 

By symmetry, the number of solutions in the open set (0, 27r) \ {7r} of the equation 



(13) 



tan — tan — = — d 
2 2 



equals the number of zeros of tan(<f^>/2) in (0, 2?r) (cf. Figure 5), which is d — 2 for even d and 
rf — 1 for odd d. Thus, all zeros of Q{d; z) are accounted for. They are either negative or complex 
conjugate pairs of zeros located on the unit circle. 

For odd d, no zero ( of G(d; z) is on the unit circle (Proposition 9) or negative (Proposition 7), 
so Q(d; C) 0, which contradicts (12). Suppose d is even. If £ is on the unit circle Co, then it is 
simple (Proposition 9). If £ is not on Co, then Q(d; () ^ 0, which again contradicts (12). □ 



*The replacements are (C - l) d (C d_1 + 1) >-» C d + C d_1 , [(C - !) d - C - l]C d_1 >-> ~(C~ l) d and £ d - (C - l) d >-» 

C d " 1 [(C-i) d -i]- 

^Interestingly, if d ^y is changed to d ^| in Q(d; z), then the new polynomial has all its zeros on the unit circle. 
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Figure 5. Typical cases when solving equation (13). 

3.3. Asymptotics of Gonchar polynomials. Numerically computing the zeros for G(d; z) for 
small values of d, we observe (cf. Figures 3 and 6 and Table 1) that they essentially form three 
groups separated by the sets 

(14a) Ai:={z£C:Rez< 1/2, \z-l\ > 1} , 

(14b) A 2 :={zeC:\z\<l,\z-l\<l}, 
(14c) A 3 :={z e C : Rez > 1/2, \z\ > 1} . 

For the purpose of asymptotic analysis (large d) we rewrite the equation G(d; z) = in three 
different ways to emphasize an exponentially decaying right-hand side when considering zeros of 
G(d; z) from the indicated part of the complex plane: 

d-1 



(15a) 
(15b) 
(15c) 



-d-1 



1 = 



1 



1 



1 



(z-l) d -(z + l)z 
(z-l) d 



d-1 



(z-l)[(z-l)z] 



d-1 



1 



(z-1) 



d-1 



(Rez < 1/2) 
-1|N<1) 
(Rez > 1/2) 



The following theorem concerning the limit behavior of the zeros of G(d; z) as d 
illustrated in Figure 6. 



oo is 



Theorem 11. Let V be the set consisting of the boundary of the union of the two unit disks 
centered at and 1 and the line-segment connecting the intersection points as indicated in Figure 6. 
Then, as d — > oo , all the zeros of G(d; z) tend to T, and every point on T attracts zeros of these 
polynomials. 



Proof. First, we observe that a closed set K in C \ T is free of zeros of G(d; z) for sufficiently large 
d as can be seen from the following relations obtained from (15a): 

l/d 



lim 

d— voo 



lim 

d— >oo 



1 



1 



Z ■ 



1 



yd-1 



1 



z- 1 



z - 1 



z-1 \z-l 

d-1 



d-1 



-d-1 



2 



lim 

d— >oo 



d-1 



l/d 



l/d 




if \z\ > 1, |z-l| > 1, 
1| if |z| > 1, \z-l\ < 1, 



if \z\ < 1 
if \z\ < 1 



|^|<l(Rez>l/2), 
|*=i| > 1 (Rez < 1/2), 



: if \z\ >1, <1 (Re* < 1/2), 

I if |z| < 1, < 1 (Rez < 1/2). 



The second part of the assertion, that every point of T attracts zeros, is proved by contradiction. 
Given a supposedly non-attracting point w on T, there is a sufficiently small 1 '' open disk D w centered 



*The disk is small enough that its intersection with T is contained either in Co, C± or the open line-segment 
connecting the intersection points of Co and Ci . 
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1 1 

-1 ~2 2 1 2 3 




-1 i i 1 2 3 



2 2 

Figure 6. Zeros of G(d; z) for d = 1, 2, . . . , 40. 

at w containing no zeros of G(d; z) for all sufficiently large d. It is possible to then define a single- 
valued analytic branch of the d-th root of any of the rational functions whose moduli appear 
on the left-hand sides above. Thereby, we obtain sequences of functions which are analytic and 
uniformly bounded in D w . Such sequences form normal families in D w . According to the right- 
hand sides above, at least one limit function of these families (which is necessarily analytic in D w ) 
will have the property that its modulus is 1 in one part and is non-constant in the other part 
of D w which is separated by T. This gives the desired contradiction, since an analytic function 
in a domain that has constant modulus on a subdomain must be constant throughout the whole 
domain. Consequently, each point of T attracts zeros of G(d; z) as d — > oo. □ 

It is inviting to compare the zeros of G(d; z) with the ones of the polynomials given at the 
left-hand sides of (15a), (15b) and (15c). Such comparisons will likely lead to a finer analysis of 
the properties of the zeros of G(d; z). 

We conclude this note with some challenging conjectures. 

Conjecture 2. For every positive integer d, the zeros of G{d; z) form three groups separated by 
the sets A\, Ai and A3 except when 6 divides d in which case one also has the zeros (1 ± iv / 3)/2. 

When counting the zeros in the sets A±, A2 and A3 a very regular pattern emerges, which can be 
seen from Table 1. In fact, inspection of this table shows that the values of column N\ (increased by 
2 when 6 divides d) agree with the number of zeros on Co obtained in Proposition 9. Assuming that 
this is true for all d > 1, by self-reprocity of G(d; 2), it would follow that Ni = N2 = N3 = 4k — 1 if 
d = 6k. We expect that the zero counting scheme indicated in the Table 1 generalizes to all d > 1. 

Numerically, the zeros in A\ and A3 can be found near the respective unit circle Cq and C\ . 

Conjecture 3. If d is even, all the zeros of G{d; z) in A\ are the zeros on the unit circle Cq 
given in Proposition 9. If d is odd, the zeros of G(d; z) in A\ alternately lie inside and outside Co- 
Furthermore, the zeros of G(d; z) in A3 are always outside of C\. 

Conjecture 4. The zeros of G{d;z) in A2 are located on a curve which is convex from the left. 
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Appendix A 

In the proof of Proposition 3 we utilize the following. 
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Lemma 12. Let d > 2 and Z\ ^ z 2 be two fixed points in R d+1 . Then 

(16) r^fi=L +1 \t^n^ dt=:J{Zl ' Z2) ' Cd = — www — • 



Proof. Observe that J(z\, z 2 ) is a convergent integral for any Zi ^ z 2 . First a translation t = u+z 2 
and then change to spherical coordinates yields 



J(zi,z 2 )= / | — \d~TTd =UJd f \ f 

J R d+l |u- Zi +z 2 | d |u| d J [J §d 



da d {u) 



dp, 



(p 2 — 2pr u • x + r 2 ) d ^ 2 j 
where we used the notation p := |u|, r := |zi — z 2 |, u := u/p, x := (zi — z 2 )/r and 

Lv d ^2^ d+1 ^ 2 /T((d+l)/2) 
is the surface area of the unit sphere § d . Using the Funk-Hecke formula [6, p. 20], we derive 



(p 2 — 2/dri; + r 2 ) c 
where the formula [7, Eq. 2.2.9.7] enables us to compute 

J ( p 2-2prv + r 2 ) d/2 d-1 V ( rf + 1 )/ 2 / 

Symmetry and a change of variable u = 1 — v 2 leads to (recall r = |zi — z 2 |) 

C d = | Z1 - z 2 |^ J( Z1 , z 2 ) = ^\[ (1 - uf 2 - 1 aPl (d - l)/2. ^ 

where the integral can be expressed as a 3F2-hypergeometric function at unity (cf. [8, Eq. 2.21.1.5]) 

- Bfd/2 1/2) ,F 2 f ^ ( d - W 2 ' d / 2 • 1 V tt [E«1±1)Z 2 11 2 

~ d—1 ' ' ' ' ' + l)/2, (<i + l)/2' y — — 1 [ r(d/2) ' 

The last step follows from the relation B(d/2, 1/2) = ujd/^d-i and the generalized hypergeometric 
function can be evaluated using [8, Eq. 7.4.4.19]. The result follows. □ 

Proof of Proposition 3. Suppose ?]i and ry 2 arc two signed equilibria on § d associated with the 
same external field Q. Then 

V m (x) + Q(x) = F 1 , V m (x) + Q(x) = F 2 for all x e § d . 

Subtracting the two equations and integrating with respect to rj := rji — r/ 2 we obtain 

l(rj) = J V^(x)drj(x) = J J | x _y |rf _ 1 dr?(x)dr ? (y) = 0. 

Applying (16) from Lemma 12 we obtain 

= / / |x-y| rf -! rf??(x) d??(y) 



d dry(x) 



1 



x-z 



dz, 



(17 » " 5 

where the interchange of the integration is justified because signed equilibria have a.e. finite 
potentials. From (17) we conclude for the Riesz potential of -q 



/-. — rj dn(x) = a.e. in 
|x-z| d 



It turns out that linear combinations of {|x — z| d } ze Rd+i are a dense class in the space of con- 
tinuous functions on S d (see [5, p. 214]), we obtain that r\ = (see also [5, Theorem 1.12]). □ 
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Proof of Asymptotic (10). Let Q denote the largest real zero of G(d,q;z) (that is R q = ( d ). We 
may assume (see discussion regarding appearance of the number 2) that 

Cd = 2 + f d with f d = f{(d) -> as d -> oo. 

Rewriting the equation G(d, q; Q) — in the following way using above relation and exploiting 
the exponential decay of the right-hand side yields 

(l + f d ) d - q (3 + f d ) = -q(l + f d )(^j^j =qO((2/3) d ) asd^co, 
or equivalently, 

log(l + f d ) = - d log(3<~) + - d log(l + ^ + G((2/3) d )) as d oo. 
Now we can use the series expansion of the logarithm function to get a relation for f d : 

fd + 0(f d ) = l - Iog(3g) + \o{f d ) as d -+ oo. 
The asymptotic form (10) follows. □ 
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